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Abstract 

We study geometric consistency relations between angles of 3-dimensional (3D) circu- 
lar quadrilateral lattices — lattices whose faces are planar quadrilaterals inscribable into 
a circle. We show that these relations generate canonical transformations of a remarkable 
"ultra-local" Poisson bracket algebra defined on discrete 2D surfaces consisting of circular 
quadrilaterals. Quantization of this structure allowed us to obtain new solutions of the 
tetrahedron equation (the 3D analog of the Yang-Baxter equation) as well as reproduce all 
those that were previously known. These solutions generate an infinite number of non-trivial 
solutions of the Yang-Baxter equation and also define integrable 3D models of statistical 
mechanics and quantum field theory. The latter can be thought of as describing quantum 
fluctuations of lattice geometry. 
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1 Introduction 



Quantum integrability is traditionally understood as a purely algebraic phenomenon. It stems 
from the Yang-Baxter equation [1,2] and other algebraic structures such as the affine quantum 
groups [3,4] (also called the quantized Kac- Moody algebras), the Virasoro algebra [5] and their 
representation theory. It is, therefore, quite interesting to learn that these algebraic structures 
also have remarkable geometric origins [6], which will be reviewed here. 

Our approach [6] is based on connections between integrable three-dimensional (3D) quantum 
systems and integrable models of 3D discrete differential geometry. The analog of the Yang- 
Baxter equation for integrable quantum systems in 3D is called the tetrahedron equation. It was 
introduced by Zamolodchikov in [7,8] (see also [9-16] for further important results in this field). 
Similarly to the Yang-Baxter equation the tetrahedron equation provides local integrability 
conditions which are not related to the size of the lattice. Therefore the same solution of the 
tetrahedron equation defines different integrable models on lattices of different size, e.g., for 
finite periodic cubic lattices. Obviously, any such three-dimensional model can be viewed as a 
two-dimensional integrable model on a square lattice, where the additional third dimension is 
treated as an internal degree of freedom. Therefore every solution of the tetrahedron equation 
provides an infinite sequence of integrable 2D models differing by the size of this "hidden third 
dimension" [9]. Then a natural question arises whether known 2D integrable models can be 
obtained in this way. Although a complete answer to this question is unknown, a few non-trivial 
examples of such correspondence have already been constructed. The first example [11] reveals 
the 3D structure of the generalized chiral Potts model [17, 18]. Another example [16] reveals 
3D structure of all two-dimensional solvable models associated with finite-dimensional highest 
weight representations for quantized affine algebras U q (sl n ), n = 2, 3, . . . , oo (where n coincides 
with the size of the hidden dimension). 

Here we unravel yet another remarkable property of the same solutions of the tetrahedron 
equation (in addition to the hidden 3D structure of the Yang-Baxter equation and quantum 
groups). We show that these solutions can be obtained from quantization of geometric inte- 
grability conditions for the 3D circular lattices — lattices whose faces are planar quadrilaterals 
inscribable into a circle. 

The 3D circular lattices were introduced [19] as a discretization of orthogonal coordinate 
systems, originating from classical works of Lame [20] and Darboux [21]. In the continuous 
case such coordinate systems are described by integrable partial differential equations (they 
are connected with the classical soliton theory [22,23]). Likewise, the quadrilateral and circular 
lattices are described by integrable difference equations. The key idea of the geometric approach 
[19,24-31] to integrability of discrete classical systems is to utilize various consistency conditions 
[32] arising from geometric relations between elements of the lattice. It is quite remarkable that 
these conditions ultimately reduce to certain incidence theorems of elementary geometry. For 
instance, the integrability conditions for the quadrilateral lattices merely reflect the fact of 
existence of the 4D Euclidean cube [25]. In Sect. 2 we present these conditions algebraically in 
a standard form of the functional tetrahedron equation [15]. The latter serves as the classical 
analog of the quantum tetrahedron equation, discussed above, and provides a connecting link 
to integrable quantum systems. 

In Sect. 13.11 we study relations between edge angles on the 3D circular quadrilateral lattices 
and show that these relations describe symplectic transformations of a remarkable "ultra-local" 
Poisson algebra on quadrilateral surfaces (see Eq. (l2ip ). Quantization of this structure allows 
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one to obtain all currently known solutions of the tetrahedron equation. They are presented in 
Sect. HI namelj0, 

(I) Zamolodchikov-Bazhanov-Baxter solution [8,11,13,34] 

(II) Bazhanov-Sergeev solution [16] 

(III) Bazhanov-Mangazeev-Sergeev solution [6] 

including their interaction-round-a-cube and vertex forms. Additional details on the corre- 
sponding solvable 3D models, in particular, on their quasi-classical limit and connections with 
geometry, can be found in [6]. 

2 Discrete differential geometry: "Existence as integrability" 

In this section we consider classical discrete integrable systems associated with the quadrilateral 
lattices. There are several ways to extract algebraic integrable systems from the geometry 
of these lattices. One approach, developed in [25,27,35-37], leads to discrete analogs of the 
Kadomtsev-Petviashvili integrable hierarchy. Here we present a different approach exploiting 
the angle geometry of the 3D quadrilateral lattices. 

2.1 Quadrilateral lattices 

Consider three-dimensional lattices, obtained by embeddings of the integer cubic lattice Z 3 into 
the iV-dimensional Euclidean space Mr, with N > 3. Let x(m) G Mr, denote coordinates of the 
lattice vertices, labeled by the 3-dimensional integer vector m = m\e\ +m2e2 +771363 G Z 3 , where 
e\ = (1,0, 0),e2 = (0,1,0) and e% = (0,0,1). Further, for any given lattice vertex xq = x(m), 
the symbols X{ = x(m + e«), Xy = x{m + Bi + ej), etc., will denote neighboring lattice vertices. 
The lattice is called quadrilateral if all its faces (xo, £Cj, Xj, Xij) are planar quadrilaterals. The 
existence of these lattices is based on the following elementary geometry fact (see Fig. 1) [25], 

3512 




Figure 1: An elementary hexahedron of a cubic quadrilateral lattice. 

Consider four points Xq, x\, X2, x% in general position in M N , N > 3. On each of 
the three planes (xQ,Xi,Xj), 1 < i < j < 3 choose an extra point not lying 

1 The other known solutions, previously found by Hietarinta [33] and Korepanov [14], were shown to be special 
cases of [34]. 



3 



on the lines (xo,Xi), (xq,xj) and (xi,Xj). Then there exist a unique point X123 
which simultaneously belongs to the three planes (xx,xi2,Xis), (X2, x\2, ^23) <ind 
(x 3 ,x n ,x 23 ). 

The six planes, referred to above, obviously lie in the same 3D subspace of the target space. 
They define a hexahedron with quadrilateral faces, shown in Fig. [TJ It has the topology of the 
cube, so we will call it "cube" , for brevity. Let us study elementary geometry relations among 
the angles of this cube. Denote the angles between the edges as in Fig. [2J Altogether we have 




Figure 2: The "front" (a) and "back" (b) faces of the cube in Fig. Q] and their angles. 



6 x 4 = 24 angles, connected by six linear relations 

OLj + Pj + 7j + Sj = 2tt, a'j + + 7J. + % = 2vr, j = 1, 2, 3, (1) 

which can be immediately solved for all "<5's". This leaves 18 angles, but only nine of them are 
independent. Indeed, a mutual arrangement (up to an overall rotation) of unit normal vectors 
to six planes in the 3D-space is determined by nine angles only. Once this arrangement is fixed 
all other angles can be calculated. Thus the nine independent angles of the three "front" faces 
of the cube, shown in Figj2^i, completely determine the angles on the three "back" faces, shown 
in Figj2}), and vice versa. So the geometry of our cube provides an invertible map for three 
triples of independent variables 

^123: {ctj,Pj,n}^{°ts,&jrfj}, j = 1,2,3. (2) 

Suppose now that all angles are known. To completely define the cube one also needs to specify 
lengths of its three edges. All the remaining edges can be then determined from simple linear 
relations. Indeed, the four sides of every quadrilateral are constrained by two relations, which 
can be conveniently presented in the matrix form 

m = *(M, x jw w)J s ^ , (3) 

V <J \ l J V C{A) D(A) ) \ isffpl »| 
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where A = {a, (3, 7,5} denotes the set of angles and £ p ,£ q ,£ p i,£ q i denote the edge lengths, 
arranged as in FigJSJ Note that due to (pQ) the entries of the two by two matrix in ([3]) satisfy 
the relation 

AD - BC = (AB -CD)/(DB - AC). (4) 
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Assume that the lengths £ p ,£ g , £ r , on one side of the two pictures in FigH] are given. Let 
us find the other three lengths £ p ',£ q ', £ r ' on their opposite side, by iterating the relation ([3]). 
Obviously, this can be done in two different ways: either using the front three faces, or the back 
ones — the results must be the same. This is exactly where the geometry gets into play. The 
results must be consistent due to the very existence of the cube in Fig. [1] as a geometric body. 
However, they will be consistent only if all geometric relations between the two sets of angles 




Figure 3: The angles A = {a, /3, 7, 6} and sides £ p ,£ q , l p < , l q ' of a quadrilateral and the oriented rapidity 

lines. 



in the front and back faces of the cube are taken into account. To write these relations in a 
convenient form we need to introduce additional notations. Note, that Figj3] shows two thin 
lines, labeled by the symbols "p" and u q" . Each line crosses a pairs of opposite edges, which we 
call "corresponding" (in the sense that they correspond to the same thin line). Eq.([3]) relates 
the lengths (£ p , £ q ) of two adjacent edges with the corresponding lengths (£' p , £' q ) on the opposite 
side of the quadrilateral. 




Figure 4: The "front" (a) and "back" (b) faces of the cube in Fig. Q] and "rapidity" lines. 



Consider now FigHk, which contains three directed thin lines connecting corresponding edges 
of the three quadrilateral faces. By the analogy with the 2D Yang-Baxter equation, where similar 
arrangements occur, we call them "rapidity" lineal. We will now apply ([3]) three times starting 
from the top face and moving against the directions of the arrows. Introduce the following three 

2 However, at the moment we do not assume any further meaning for these lines apart from using them as a 
convenient way of labeling to the corresponding (opposite) edges of quadrilaterals. 
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by three matrices 

/ A B \ / A B \ / 1 \ 

X p? (.A) = C £> , X pr (.4) = 1 , X qr (A) = [ A B \ , (5) 

\ 1 / \ C D J \ C D J 

where A,B,C,D are defined in ([3]) and their dependence on the angles A = {a, /3, 7, 5} is 
implicitly understood. It follows that 

(^.^.O* = X M (A)^pr(^2)^(^3) (ip^irf (6) 

where 

A' = {aj>#7>7j><*j}) i = 1 , 2 ,3, (7) 

the lengths i p ,£„,... are defined as in FigU and the superscript "t" denotes the matrix trans- 
position. Performing similar calculations for the back faces in FigHJa and equating the resulting 
three by three matrices, one obtains 

X pq {Ai) X pr {A2) X qr {Az) = X qr (A' 3 ) X pr (A 2 ) X pq (A[) . (8) 

where 

4 = {c4,$, 7 ^}, j = 1,2,3. (9) 

This matrix relation contains exactly nine scalar equations where the LHS only depends on the 
front angles ([7]), while the RHS only depends on the back angles ([9]). Solving these equations 
one can obtain explicit form of the map ([2]). The resulting expressions are rather complicated 
and not particularly useful. However the mere fact that the map ([2]) satisfy a very special Eq.([8]) 
is extremely important. Indeed, rewrite this equation as 

X pq (Ai) X pr (A 2 ) X qr (A 3 ) = K 123 (XgriAs) X pr (A 2 ) X pq (Aifj (10) 

where TZ\23 is an operator acting as the substitution ([2]) for any function F(Ai,A2,As) of the 
angles, 

Then, following the arguments of [14], one can show that the map ([2]) satisfies the functional 
tetrahedron equation [15] 

7^-123 ' 7^145 • 7^-246 ' ^.356 = ^-356 • 7^246 ' ^145 • 7^-123 , (12) 

where both sides are compositions of the maps ([2]), involving six different sets of angles. Alge- 
braically, this equation arises as an associativity condition for the cubic algebra (fTU|) . To discuss 
its geometric meaning we need to introduce discrete evolution systems associated with the map 
©• 

2.2 Discrete evolution systems: "Existence as integrability" 

Consider a sub-lattice L of the 3D quadrilateral lattice, which only includes points x{m) with 
mi, ni2, m?, > 0. The boundary of this sub-lattice is a 2D discrete surface formed by quadrilat- 
erals with the vertices x(m) having at least one of their integer coordinates mi, m 2 , m 3 equal to 
zero and the other two non-negative. Assume that all quadrilateral angles on this surface are 
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known, and consider them as initial data. Then repeatedly applying the map one can calcu- 
late angles on all faces of the sub-lattice L, defined above (one has to start from the corner x(0)). 
The process can be visualized as an evolution of the initial data surface where every transforma- 
tion (|2|) corresponds to a "flip" between the front and back faces (Fig. [2]) of some cube adjacent 
to the surface. This makes the surface looking as a 3D "staircase" (or a pile of cubes) in the 
intersection corner of the three coordinate planes, see Fig. [5] showing two stages of this process. 
Note, that the corresponding evolution equations can be written in a covariant form for an arbi- 



trary lattice cube (see Eq. (|23p below for an example). It is also useful to have in mind that the 
above evolution can be defined purely geometrically as a ruler- and- compass type construction. 
Indeed the construction of the point #123 in Fig. [T]from the points Xq, x%, x%, X3, x\2, £Ci3, X23 
(and that is what is necessary for flipping a cube) only requires a 2D-ruler which allows to draw 
planes through any three non-collinear points in the Euclidean space. 

Similar evolution systems can be defined for other quadrilateral lattices instead of the 3D 
cubic lattice considered above. Since the evolution is local (only one cube is flipped at a time) 
one could consider finite lattices as well. For example ,consider six adjacent quadrilateral faces 
covering the front surface of the rhombic dodecahedron^ shown in Fig. [6j Suppose that all angles 
on these faces are given and consider them as initial data. Now apply a sequence of four maps 
(|2|) and calculate angles on the back surface of the rhombic dodecahedron. This can be done in 
two alternative ways, corresponding to the two different dissection of the rhombic dodecahedron 
into four cubes shown in FigJBJ The functional tetrahedron equation (|12p states that the results 
will be the same. Thereby it gives an algebraic proof for the equivalence of two "ruler-and- 
compass" type constructions of the back surface of the dodecahedron in Fig. [6J Can we also 
prove this equivalence geometrically? Although from the first sight this does not look trivial, 
it could be easily done from the point of view of the 4D geometry. The required statement 
follows just from the fact of existence of the quadrilateral lattice with the topology of the 4D 
cube [25]. The latter is defined by eight intersecting 3-planes in a general position in the 4-space. 
The two rhombic dodecahedra shown in Fig. [6] are obtained by a dissection of the 3-surface of 
this 4-cube, along its 2-faces, so these dodecahedra must have exactly the same quadrilateral 
2-surface. Thus the functional tetrahedron equation (|12p . which plays the role of integrability 
condition for the discrete evolution system associated with the map ([2]), simply follows from the 
mere fact of existence of the 4-cube, which is the simplest 4D quadrilateral lattice. For a further 

3 It is worth noting that the most general rhombic dodecahedron with quadrilateral faces can only be embedded 
into (at least) the 4D Euclidean space. 




Figure 5: Visualization of the 3D "staircase" evolution. 
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Figure 6: Two dissections of the rhombic dodecahedron into four quadrilateral hexahedra. 



discussion of a relationship between the geometric consistency and integrability see [32]. 

3 Quantization of the 3D circular lattices 
3.1 Poisson structure of circular lattices 

The 3D circular lattice [19,28,29] is a special 3D quadrilateral lattice where all faces are circular 
quadrilaterals (i.e., quadrilaterals which can be inscribed into a circle). The existence of these 
lattices is established by the following beautiful geometry theorem due to Miquel [38] (see Fig. [7]) 




Figure 7: Miquel configuration of circles in 3D space, an elementary hexahedron and its circumsphere. 
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Miquel theorem. Consider four points Xq, X\, X2, X3 in general position in 
N > 3. On each of the three circles c(xo,Xi,Xj), 1 < i < j < 3 choose an additional 
new point Xij. Then there exist a unique point 2:123 which simultaneously belongs to 
the three circles c(a?i, X12, X13), c(£E2, 2:12, #23) and c{x^, #13, X23). 

It is easy to see that the above six circles lie on the same sphere. It follows then that every 
elementary "cube" on a circular lattice (whose vertices are at the circle intersection points) is 
inscribable into a sphere, see Fig. [7J The general formulae of the previous subsection can be 
readily specialized for the circular lattices. A circular quadrilateral has only two independent 
angles. In the notation of Fig. [3] one has 



7 = 7T - P, 



it — a . 



(13) 



Due to the Miquel theorem we can simply impose these restrictions on all faces of the lattice 
without running to any contradictions. The two by two matrix in ([3|) takes the form 



X 



\ -a k J 
where we have introduced new variables 

k = (esc a) sin f3 , a = (esc a) sin(a + 0) 



det X = 1 



(esc a) sin(a — /?), 



(14) 



(15) 



instead of the two angles {a, ft}. Note that the new variables are constrained by the relation 



Conversely, one has 



aa 



a — a 



cos a 



2k 



1-k 2 



cos (3 



a + a* 



(16) 



(17) 



Let the variables {kj, dj, a*}, {k'j, a'-, af}, j = 1,2,3, correspond to the front and back faces of 
the cube. The map ([2D then read explicitly 



n 



123 



where e = +1 and 



^3^1 ^ 1 ^2^3 ' 



{k 2 a\)' 
(a 2 )' = a\a\ + ek^a^ 



(^2 a l)' — ^3 a l e ^l a 2 a 3> 

(a 2 )' = 01^3 + ek 1 k 3 a 2 , 



(18) 



( ^3 ) — ^1^3 sk^d^ CLq ) (^2^3) — ^1^3 sk^Cl^Qi- 



2- 



1 - a' 2 a* 2 ' 



(19) 



At this point we note that exactly the same map together with the corresponding equations 
([8]) and ([T2j) were previously obtained in [16]. Moreover, it was discovered that this map is a 
canonical transformation preserving the Poisson algebra 



} — 2 5{j fcj , {ki,aj} — 5ij ki aj , {fej,a„-} 



5ij ki dj , 



hi = 1,2,3 



(20) 



where kf = 1 — OjO*. Note that variables k,a,a* on different quadrilaterals are in involution. 
The same Poisson algebra in terms of angle variables reads 



{ai,(3j} = Sy , {a i7 aj} = = 



(21) 
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This "ultra-local" symplectic structure trivially extends to any circular quad-surface of initial 
data, discussed above. To resolve an apparent ambiguity in naming of the angles, this surface 
must be equipped with oriented rapidity lines, similar to those in Fig. In addition, the angles 
for each quadrilateral should be arranged as in Fig. [3l Then one can assume that the indices 
i, j in (12ip refer to all quadrilaterals on this surface. 

Thus, the evolution defined by the map (|18p is a symplectic transformation. The correspond- 
ing equations of motion for the whole lattice (the analog of the Hamilton-Jacobi equations) can 
be written in a "covariant" form. For every cube define 

^32 = a l> ^23 = a l; ^31 = a 2; ^13 = a 2*> ^21 = a 3; ^12 = a 3 (22) 

where Ajj. stands for Ajk(m), where m is such that x{m) coincides with the coordinates of the 
top front corner of the cube (vertex xq in FigfT]). Let T k be the shift operator TkAij(m) = 
Aij(m + ek)- Then 

An AjuA 



TkA{j — , K%j — Kji — yl AijAji , (23) 

where (i,j,k) is an arbitrary permutation of (1,2,3) and 

f x =T x , f 2 = T 2 -\ T 3 = T 3 . (24) 

Note that Eq. (f2"3"J) also imply 

{T k K i:j )K kj = {TiK^Kij . (25) 

Remarks. The equations (|23p have been previously obtained in [29], see Eq.(7.20) therein. 
The quantities A{j in (|23p should be identified with the rotation coefficients denoted as f5ij 
in [29]. The same equations (|23p are discussed in §3.1 of [32], where one can also find a de- 
tailed bibliography on the circular lattices (we are indebted to A.I.Bobenko for these important 
remarks) . 



3.2 Quantization and the tetrahedron equation 

In the next section we consider different quantizations of the map (|18p and obtain several 
solutions of the full quantum tetrahedron equation (see Eq. (|33p below). In all cases we start 
with the canonical quantization of the Poisson algebra (|2ip . 

[ai,Pj] = [a l ,a j }=0, [f3 h fa] = , (26) 

where h is the quantum parameter (the Planck constant) and £ is a numerical coefficient, in- 
troduced for a further convenience. The indices i,j label the faces of the "surface of initial 
data" discussed above. Since the commutation relations (|26p are ultra-local (in the sense that 
the angle variables on different faces commute with each other), let us concentrate on the local 
Heisenberg algebra, 

H: [a,0\=th, (27) 

for a single lattice face (remind that the angles shown in Fig. [3] are related by (|13p ). The 
map (fT8l) contains the quantities k,a,a*, defined in (fToT) . which now become operators. For 

4 We refer the reader to our previous paper [39] where the relationship between the rapidity graphs and 
quadrilateral lattices is thoroughly discussed. 
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definiteness, assume that the non-commuting factors in (|15|) are ordered exactly as written. 
Then the definitions (|15p give 



k = (U - U- 1 )- 1 (V - V- 1 ), 

a = q~h (U - I!' 1 )- 1 (UV — U' 1 V^ 1 ), (28) 

a* = q+\ (U - If' 1 )- 1 {U V~ l - U- 1 V), 
where the elements U and V generate the Weyl algebra, 

UV = qVU, U = e ia , V = e if3 , q = e~ ih . (29) 
The operators (|28p obey the commutation relations of the (/-oscillator algebra, 



OsCq : 



qa*a — q 1 aa*=q — q , ka*=qa*k, ka = q 1 ak 
k 2 = q (1 - a* a) = q' 1 (1 - a a*) . 



(30) 



where the element k is assumed to be invertible. This algebra is, obviously, a quantum counter- 
part of the Poisson algebra (|20p . In the previous Section we have already mentioned the result 
of [16] that 

(i) the map (|18p is an automorphism of the tensor cube of the Poisson algebra (|20p (remind 
that the relation (TTBI) should be taken into account in (fT8l)). 



In the same paper [16] it was also shown that 

(ii) there exists a quantum version of the map (|18|) . which acts as an automorphism of the 
tensor cube of the g-oscillator algebra (|30p . The formulae (|18p for the quantum map stay 
exactly the same, but the relation (|16p should be replaced by either of the two relations 
on the second line of ([30]) . for instance, k 2 = q (1 — a* a). In particular, (fT9|) should be 
replaced with 

(k' 2 ) 2 = q (l-a* 2 'a 2 ') . (31) 

(hi) the quantum version of the map (|18p . defined in (ii) above, satisfies the functional tetra- 
hedron equation (fT2|) . 



For any irreducible representations of the g-oscillator algebra (|30|) the formulae (|T8|) and (|3T|) 
uniquely determine the 7^123 as an internal automorphism, 

^123 {F) = R 123 F R u \, F G Os C(? 8> Oscg ® Osc q . (32) 

It follows then from (|12p that the linear operator -R123 satisfies the quantum tetrahedron equation 

-Rl23 -Rl45 -^246 -^356 = -^356 -^246 ^145 -^123 j (33) 

where each of the operators -R123, R145, R246 and -R356 act as (|32p in the three factors (indicated 
by the subscripts) of a tensor product of six g-oscillator algebras and act as the unit operator 
in the remaining three factors. 
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4 Solutions of the tetrahedron equation 



Here we show that all known solutions of the tetrahedron equation can be obtained by solving 
(|32p for the linear operator -R123 with different irreducible representations of the q-oscillator 
algebra (1501) . 

4.1 Fock representation solution 

In this subsection we set £ = 1 in (|26p and q = e~ h . Define the Fock representation of a single 



g-oscillator algebra (|30|) . 

a|0)=0, a\n+l) = \n), a* \n) = (1 - q 2+2n ) \n + 1), k \n) = q n+l / 2 \ n ), (34) 

spanned by the vectors |n), n = 0, 1, . . . , 00. Then using ([18]) . ([30]) . (f5T|) and ([32|) one can show 
that 



(711,712, n 3 |fl|ni, n' 2 ,n 3 > = <5„ 1 +„ 2 ,„ i +„^„ 2 +„ 3 ,„ i +^ (-1)" 2 q {n{ nM ™ a) 



x H) 2Mq- 2n ' 2 ,q 2{1+n '"\q 2{1 - n ' 2+n3) -,Q 2 ,Q 2{1+ni) ) , 



(35) 



where 



and 



-nc -«"»>• fl) -rafc. < 36 > 



i=o 



, / , 2 \ («;g 2 )«(M 2 )n „ 

*Ma,b,c;q ,Z) = ^ iq2]q2Uc;q2)n Z (37) 

is the q-deformed Gauss hypergeometric series. This 3D f2-matrix satisfies the constant tetra- 
hedron equation (|33p . In matrix form this equation reads 



Ejjn 1 ,n 2 ,n 3 p n l ' n 4' ,i 5 p n 2> n 4> n 6 p n 3' n 5> n 6 

n'.=0 
3 

oo 

Er> n 3' n 5< n 6 r> n 2< n 4> n 6 rXlOs R n i>2>3 
-rtn 3 ,n 5 ,n 6 ^^^J, 

nj=0 



(38) 



where the sum is taken over six indices n^, n' 2 , n' 3 , n'^, n' 5 , n 6 and 

= (ni,n2,n 3 \R\n' 1 ,n 2 ,n' 3 ) . (39) 

Note that Eq. (|38p does not contain any spectral parameters. Originally the i?-matrix (|35p was 
obtained in [16] in terms of a solution of some recurrence relation, which was subsequently 
reduced to the g-hypergemetric function in [6]. 

4.2 Modular double solution 

In this subsection we set in (|26p 

£ = -t, h = 7rb 2 , q = e i7Tb2 (40) 
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where b is a freeparameter, Re b ^ 0. Here it will be more convenient to work with a slightly 
modified versioqfl of the map (|18p . with the value e = — 1. 

Consider a non-compact representation [40] of the g-oscillator algebra (|30p in the space of 
functions /(<r) G L 2 (IR) on the real line admitting an analytical continuation into an appropriate 
horizontal strip, containing the real axis in the complex cr-plane (see [40] for further details). 
Such representation essentially reduces to that of the Weyl algebra 

W 9 : kw = qwk, q = e™ h \ (41) 

realized as multiplication and shift operators 

k \a) = ie 7 ™" 6 \a) , w\a) = \a-\b). (42) 

The generators a, a* in (f30|) are expressed as 

a = w~ 1 , a* = {l-q- 1 k 2 )w . (43) 

As explained in [41] the representation (I42j) is not, in general, irreducible. Therefore, the relation 
(|32p alone does not unambiguously define the linear operator -R123 in this case. Following the 
idea of [41] consider the modular dual of the algebra ([4"Ijh 

W g ~: kw = qwk, q = e - i7Tb ~ 2 , (44) 

acting in the same representation space 

k \a) = -ie™ 6 " 1 \a) , w \a) = \a + i6 _1 ). (45) 

We found that if the relation (|32p is complemented by its modular dual 

£123 [F) = R123 F Ru3, F e Osc^ ® Oscg ® Osc 9 - , (46) 

then the pair of relations (f32|) and (j46]l determine the operator -R123 uniquehj^. The dual (7- 



oscillator algebra Oscq is realized through the dual Weyl pair (|44p and the relations 

5 = u}" 1 , a* = (1 - g-ip) t5 . (47) 

The dual version of the map 7^123 is defined by the same formulae (fl~8|) . where quantities kj, aj,a*j, 
j = 1, 2, 3 are replaced by their "tilded" counterparts kj,aj, a*. The value of q does not, actually, 
enter the map (fTHf) . but needs to be taken into account in the relations between the generators 
of the g-oscillator algebra. Thus, the linear operator -R123 in this case simultaneously provides 
the two maps 7^123 and 7^123 (given by (fT8l) with e = — 1). The explicit form of this operator is 
given below. 
Denote 



and define a special function 

a+!ZZ) y(z + £2+12) 



2^2 



f ci '%) = / ^ e ^(- C0 -i,)^±zizMi+-T- 

V c 3,c 4 J Jr l n( z + %JR) lp (z + 



(49) 



5 Note that Eq. (|18[) is a particular case a general map considered in [16]; see also Eq. (|59|l below. 

6 It is worth mentioning similar phenomena in the construction of the i?-matrix [42] for the modular double of 



the quantum group U q (sl2) and the representation theory of U q (sl2,M.) [43]. 
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where (p is the non-compact quantum dilogarithm [44] 



<p(z) 



ex Pl 4 



-2\zx 



dx\ 

, i0 sinh(j;6)sinh(x/6) x J 



(50) 



The values of ci, 02,03,04 are assumed to be such that poles of numerator in the integrand of 
(|49p lie above the real axis, while the zeroes of the denominator lie below the real axis. For other 
values of Cj the integral (|49[) is defined by an analytic continuation. Note that for Im b 2 > the 
integral 2^2 can be evaluated by closing the integration contour in the upper half plane (see 
eq.(78) in [6]), which is very convenient for numerical calculations. 

The matrix elements of i?-matrix solving the pair of the equations (|32p and (|46p are given 

by 



(01, a 2 ,a 3 \R\a[,a' 2 ,a' 3 ) = 5 a 



T> X 



x e 



i7r K "3+ ir ?(°"l+ cr 3- cr 2)) 



2^2 



01—03, —01 + 03 

01 + 03, -<A - 3 



02 



(51) 



This ii-matrix satisfies the constant tetrahedron equation f|33|) . Its matrix form is given by (|38j) 

where i4<!n^nfc is substituted by i?<T-^'o- fc = (0i, 0j, 0fe|-R |cr^, ct^-, cr^.), i, j, = 1, 2, . . . , 6, and the 
sums are replaced by the integrals over a[, a' 2 , cr' 3 , a' 4 , a' 5 , a' 6 along the real lines — 00 < a\ < 00. 
One can verify that these integrals converge. The solution (|5ip was obtained in [6]. 



4.2.1 The "interaction-round-a-cube" formulation of the modular double solution 



Note that due to the presence of two delta-functions in (I5ip the edge spins are constrained by 
two relations a 1 + a 3 = a[ + a 3 and a 2 + 3 = 2 + 03 at each vertex of the lattice. Here we 
re-formulate this solution in terms of unconstrained corner spins, which also take continuous 
values on the real line. Figure [8] shows an elementary cube of the lattice with the corner spins 
li a,b,c,d,e, f, g,h" arranged in the same way as in [10]. The corresponding Boltzmann weight 
reads [6] 




Figure 8: The arrangement of corner spins around a cube. 
W(a\e, f,g\b,c,d\h;Ji,T 2 ,T 3 ) = 

>«4 + ir,K+4-<x 2 )) ^ 2 f 01 " 03, -01 + 03 

01+03, —0i - 03 



02 



where 



01 = 9 + / — a — b — Ti, a 2 = a + c 



o + T 2 







c + d 



h-Tt, a' 2 = f + h-b- d + J 



2- 



03 



(To 



e + /- a- d-T 3 , 
b + c- g -h-J 3 



(52) 



(53) 
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Note, that the variables <7j and a\ automatically satisfy the delta function constrains in (|5ip . 
The weight functions (|52p contains three arbitrary (spectral) parameters Ti,T2,T3. It satisfies 




Figure 9: Graphical representation for the tetrahedron equations for interaction-round- a-cube models, 
the tetrahedron equation in the "interaction-round-a-cube" form [10] (see Fig. [9j), 
dz W(a 4 \d,c 3 , c 2 \b 3 , b 2 , h\z) W'(cx\a 3 , h,b 2 \z, c 6 , c 4 |6 4 ) 
xW"(bi\c4, c 3 , z\b 3 , 64, a 2 |c 5 ) W"'(z\h, b 3 , b 2 \c 2 , c 6 , c 5 \ai) = 



J 



(54) 



dz W'"(bi\c4, c 3 , ci, |a 4 , a 3 , a 2 |z) ^"(^103,04, b 2 \c 2 , c 6 , z|ai) 

x ^'(04^,03,02163,01,02105) VF(z|a 3 ,a2,ai|c5,c 6 ,C4|64), 

where the four sets of the spectral parameters are constrained as 

T i = T i> T 1 = -T 2 , T x =T 3 , T 2 =T 2 , T 2 = -T 3 , T 3 = T 3 . (55) 

Note that the parameters Ti,T2,T3 are similar to those in the Zamolodchikov model [8] and its 
generalization for an arbitrary number of spin states [11], which is considered in Sect. 14.31 below. 
One can relate them as Tj = log[tan(#j/2)] to the angles 9i,9 2 ,6 3 of the spherical triangle in 
Sect. 14.3721 In total, there are twelve parameters Jj, T'- ; "Pj, T'j', j = 1, 2, 3 (three for each weight 
function) constrained by six relations (|55p . Thus, the tetrahedron equation (|54p contains six 
independent parameters (in contract to only five parameters in the cyclic case; see the text after 
Eq.([72D below). 

The above solution can be simply generalized by multiplying the weight (I52|) by an "external 
field" factor 

3 

W(. . .) -» exp fjivj + ty] W(. . .) . (56) 
i=i 

This substitution does not affect the validity of the tetrahedron equation (|54p , provided the field 
parameters fj for different W's are constrained as 

rll rl r rill rll rl rill rll , r rill r rl /K^7\ 

J3—J3-J31 h =h~ /li h — h + h, J 3 —72-/2- 

Similar generalizations apply for the solutions (|35|) and (|5ip . The solution (|52p was previously 
obtained in [6]. 
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4.3 Cyclic representation solution 

4.3.1 Generalized form of the (/-oscillator map 



Let 

\ 





with a ] 


( i 








Xjkj 







V o 






aj 



J = 1,2,3, (58) 



-Hjkj 
Xjfij ) 

be an operator-valued matrix acting in the direct product of two vector spaces C 2 ® C 2 , whose 
elements depends on a set of generators Hj = {kj,aj,a*} of the (/-oscillator algebra (f30l) and 
continuous parameters A,-, and fij, where j = 1,2,3. The RHS of (|58|) is understood as a block 
matrix with two-dimensional blocks where matrix indices of the second space C 2 numerate the 
blocks, while those for first space numerate the elements inside the blocks. 
The required generalization of (|32[) reads [16jj| 

Ln(Hi) L 13 (H 2 )L 23 (H 3 ) = R123 ^23(^3) L 13 (H 2 ) L l2 {H x ) R u \ (59) 



Here £12 denote the matrix which acts as (158j) in the first and second component of the tensor 
product C 2 ® C 2 ® C 2 and coincides with the unit operator in the third component (£13 and £23 
are defined similarly). It is important to note that the additional parameters Aj and pj enter 
(I59|) only in three combinations A2/A3, X±p 3 and pi/p 2 (explicit formulae for the corresponding 
generalization of the map (|18p are given in Eqs. (24-26) of [16]). 

4.3.2 Cyclic representations of the g-oscillator algebra 

Let N > 3 be an odd integer, 

q = - e ™/ N , q N = 1, N = odd, N > 3. (60) 

and V N be an TV-dimensional vector space, spanned by the vectors |n), n G Z^r = {0, 1, . . . , N — 
1}. Define N by N matrices (the index n is treated modulo N) 

X\n)=q n \n), Z\n) = \n + 1), X N = Z N = 1, (61) 

and their embedding into a direct product V N (g) V N ® V , 

Xj = l®---® X = 1® •••® Z j = 1,2,3. (62) 

i-th i-th 



Equation (|59j) involves the direct product of three copies of the g-oscillator algebra (l30j) . labeled 
by the subscript j = 1,2,3. The most general cyclic representation for this product 

k^XjXj, a* = -(\-q- l x)X])Z 3 , a, = Pj Zj\ j= 1,2,3, (63) 
Pi 

contains two continuous parameters Xj,Pj, for each factor, so there six arbitrary parameters 
altogether. Detailed inspection of (|59|) shows, however, that pi,p 2 ,p 3 only enter through a ratio 
P1P3/P2, so there are only four essential parameters. 



Equation (I59|l reduces to 11321) when Ai = A2 = \s = 1 and jiti = /U2 = H3 = —£, where e enters the map (|18|) 
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Let 61, #2; #3 and 01,02,03 be angles and sides of a spherical triangle, respectively. Define 
2/3i = a 2 + a 3 - a%, 2(3 2 = ai + a 3 - a 2 , 2/3 3 = ai + a 2 - a 3 , /? = vr - /3i - /3 2 - /3 3 . (64) 
and set 

xi = i ^/tan y, x 2 = i yjcot y, x 3 = i ^/tany . (65) 
Then Eq, (|59p implies 

P1P3/P2 = e~ ] ^ N Vsmo 2 /sin/3 , (66) 

and 

A 3 /A 2 = e 1 ^ , A 1M3 = e-W , fi 1 /p 2 = e' a ^ N . (67) 

4.3.3 Solution of the tetrahedron equation 

Let p denote a point on Fermat curve 

p = (x,y I x N + y N = 1). (68) 
Introduce a function <p p (n), n £ Z^r, (cyclic analog of the quantum dilogarithm) [11] 

^(0) = 1, ^ ( "~ . 1} = l^gg^ ^(n + iV) = ^(n) (69) 

The i?-matrix solving (|59l) for the cyclic representation (|63p is given by 

{n\R\n } = 5 ni+n2 ^ n ' i+n ^5 n , 2+n3 ^ +n i ri x 

funa-nJ(ni-H»s) g -2nn^ ^Pi ( n + n l + n 3)^P2 M (70) 

7^ Vp 3 ( ra + "-1)^4(^ + ^3) ' 

where the points pj on the curve (|68p are defined by 

xi = e"^ Vsm/3 2 /sin/3 , y x = e'^ N Vsina 2 /sin/% , 

X2 = e -ia 2 /iv y sin/3o / sin/ 3 2 , y2 = e ift>/* V s ina 2 /sin/?2 , 

x 3 = e- 1 ^ 2 ^)/^ Vsin^s/sin/?! , y 3 = e~^ N Vsina 2 /sin/?i , 
X4 = e - i ( a 2+^)/^v ^sin^/sin^ , y 4 = e"^ 1 ^ ysina 2 /sin/3 3 . 



(71) 



We write this i?-matrix as i?i23(^i, 02 > 63) to indicate its dependence on the angles 01,02,03- It 
satisfies the tetrahedron equation (133p . provided the i?-matrices therein are parametrized by 
four triples of dihedral angles at four vertices of an Euclidean tetrahedron 

-Rl23 = ^123(01)02) #3)5 -^145 = ^145(01) 7T — 04, 7T — 0s), 

(72) 

-^246 — -R246( 7r — #2) 7T — 6*4, 0q), -^356 — -^356(03) 05) 04) • 

Here 0, , j = 1, 2, . . . , 6 are inner dihedral angles of the tetrahedron. These angles are constrained 
by one relation, therefore Eq. (|33p in this case contains five arbitrary parameters. Finally, the 
substitution 

n\ = g + f — a — 6, n 2 = a + c — e — g n 3 = e + / — a — d, 

(73) 

n\ = c + d — e — h, n' 2 = f + h — b — d, n 3 = 6 + c — g — h 
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and some trivial equivalence transformations bring (|7U|) into its "interaction-round-a-cube" form 
W(a\eJ, 9 \b,c, d \h) = £ g 2n( W -,) ^(n-^cH 2 (n-/ + ^ 

which satisfies (|54p where the integration is replaced by the summation over Z^v, while the angle 
parametrization remains the same as in (|72|) . Although our considerations in the cyclic case 
were restricted to odd values of N in (|60p . the final result (|74p as valid for even values of N as 
well. The solution (|74j) was previously obtained in [11]; for N = 2 it coincides with the Baxter's 
form [10] for the Boltzmann weights the Zamolodchikov model [8]. 



5 Conclusion 

In this paper we have exposed various connections between discrete differential geometry and 
statistical mechanics, displaying geometric origins of algebraic structures underlying integrability 
of quantum systems. 

We have shown that the 3D circular lattices are associated with an integrable discrete Hamil- 
tonian system and constructed three different quantizations of this system. In this way in Sect.H] 
we have obtained all previously known solutions of the tetrahedron equation. They are given by 
Eqs. (|35|) . (j5Tj) . ([52]) . (|70l) . The resulting 3D integrable models can be thought of as describing 
quantum fluctuations of the lattice geometry. The classical geometry of the 3D circular lattices 
reveals itself [6] as a stationary configuration giving the leading contribution to the partition 
function of the quantum system in the quasi-classical limit. 

The solutions of the tetrahedron equation discussed here possess a remarkable property: 
they can be used to obtain [11,16] infinite number of two-dimensional solvable models related to 
various representations of quantized affine algebras U q (sl n ), n = 2, 3, . . . , oo (where n coincides 
with the size of the "hidden third dimension"). Apparently, a similar 3D interpretation, origi- 
nating from other simple geometrical models, also exists for the trigonometric solutions of the 
Yang-Baxter equation, related with all other infinite series of quantized affine algebras [45,46] 
and super-algebras [47] (see [48,49] for recent results in this direction). Therefore, it might 
very well be that not only the phenomenon of quantum integrability but the quantized algebras 
themselves are deeply connected with geometry. 

There are many important outstanding questions remain, in particular, the geometric mean- 
ing of the Poisson algebra (|2ip and connections of the 3D circular lattices with the 2D circle 
patterns [50] on the plane or the sphere. It would also be interesting to understand underlying 
reasons of a "persistent" appearance of the g-oscillator algebra ([30]) as a primary algebraic struc- 
ture in many other important aspects of the theory of integrable systems, such as, for example, 
the construction of Baxter's Q-operators [51] and the calculation of correlation functions of the 
XXZ model [52]. It would also be interesting to explore connections of our results with the in- 
variants of the 3D manifolds [53-55], the link invariants [56-58], quantization of the Techmueller 
space [59,60] and the representation theory of U q (sl(2\M.) [61]. So, there are many interesting 
questions about the quantum integrability still remain unanswered, but one thing is getting is 
more and more clear: it is not just connected with geometry, it is geometry itself! (though the 
Quantum Geometry). 
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